ac Josephson effect in finite-length nanowire junctions with Majorana modes 
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It has been predicted that superconducting junctions made with topological nanowires hosting 
Majorana bound states (MBS) exhibit an anomalous 47r-periodic Josephson effect. Finding an 
experimental setup with these unconventional properties poses, however, a serious challenge: for 
finite-length wires, the equilibrium supercurrents are always 27r-periodic as anticrossings of states 
with the same fermionic parity are possible. We show, however, that the anomaly survives in 
the transient regime of the ac Josephson effect. Transients are moreover protected against decay 
by quasiparticle poisoning as a consequence of the quantum Zeno effect, which fixes the parity of 
Majorana qubits. The resulting long-lived ac Josephson transients may be effectively used to detect 
MBS. 

PACS numbers: 74.50. +r 73.21.Hb 03.65.Vf 



Recently, it has been argued that Majorana bound 
states (MBS) should appear in topological insulators [1 
and semiconductors with strong spin-orbit (SO) cou- 
pling 2-5 which, in proximity to s-wave superconduc- 
tors, may behave as topological superconductors (TS). 
MBS in these TS can be understood as Bogoliubov-De 
Gennes (BdG) quasiparticles appearing inside the super- 
conducting gap, exactly at zero energy (for a review see 
Ref. [5]). The TS phase is tunable, which has spurred 
a great deal of experimental activity towards detecting 
MBS. Tunneling through such TS is expected to show 
a zero-bias anomaly signaling the presence of MBS [TJ- 
113] . However, such anomaly, which has been recently 
observed [T3] . only proves a necessary condition. A 
a more stringent test can be established by measuring 
the Josephson effect through a junction between two TS 
nanowires. Kitaev predicted [15] that such Josephson 
effect has an anomalous 4-7r-periodicity in the supercon- 
ducting (SC) phase difference 4> = <t>i ~ §i between the 
two wires. This fractional Josephson current Ij ~ sin(^) 
[4ll5l fT5H20] can be understood in terms of fermion parity 
(FP): if FP is preserved, there is a protected crossing of 
Majorana states at <j> = tt [21] with perfect population in- 
version, namely the system cannot remain in the ground 
state as cf> evolves from to 2ir adiabatically, unlike for 
standard Andreev bound states (ABS) [52]. However, 
any finite length of the two TS regions gives rise to two 
additional MBS at the wire ends (which are connected to 
topologically trivial regions), allowing for the hybridiza- 
tion of two states of the same FP. This results in residual 
splittings at <p = it which, despite being exponentially 
small, destroy the fractional effect as the system remains 
in the ground state for all </> [23] . 

The 47r-periodicity can be restored in two ways. On 
the one hand, one may employ a fully topological circuit 
[2~4"] . which is free of additional MBS. In other words, by 
employing an S'NS' geometry (where S' is a TS and N 
is a normal, non-topological region), as opposed to an 
SS'NS'S junction (where S is a topologically trivial su- 
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FIG. 1. (Color online) (a) Schematics: A nanowire of length 
L = Lgi + Ln + L$> in contact with two s-wave supercon- 
ductors (with gap As), develops a proximity-induced SC gap 
As' < As and four Majorana modes 71,2,3,4. The trans- 
parency of the junction can be controlled by a gate voltage 
V g . (b) Energy of the four lowest many-body states in the 
TS phase, as a function of SC phase difference 4> across A?". 
They correspond to even (solid) and odd (dashed) fillings of 
the two lowest ABS, see Fig. [2]3. Note the avoided crossing 
of size St, at (j> — n, and the detachment gap So at (f> — 0. 



perconductor). Unfortunately, this is presently a difficult 
experimental challenge using nanowires. Alternatively, 
the 47r Josephson effect may be recovered by biasing the 
junction and thereby sweeping <fi fast enough through the 
anticrossing such that Landau-Zener (LZ) processes [25] 
induce non-adiabatic transitions between states of the 
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same FP. By the same token, however, transitions to the 
continuum of states above the gap, and hence parity mix- 
ing processes, will be unavoidable [26] . 

The above hindrances pose relevant questions. Is the 
Josephson effect a useful tool to detect MBS in realistic 
TS wires? Do finite TS length and quasiparticle poi- 
soning inevitably destroy the fractional periodicity? We 
here prove that the ac currents in the TS phase con- 
tain anomalously long-lived 47r-periodic transients which 
are tunable through both bias and gate voltages. Thus, 
the Josephson effect may still be exploited to provide 
unequivocal proof of the existence of MBS in finite- 
length TS. Interestingly, the duration of the long 4-7T- 
periodic transients may be increased by the well-known, 
though somewhat counter-intuitive, quantum Zeno effect, 
whereby a strong coupling to a decohering environment 
helps to confine the dynamics of a quantum system into 
a desired sector of the Hilbert space [27] . 

Low-energy description.- Consider a semiconducting 
one-dimensional wire with chemical potential fig> i SO 
coupling a, and Zeeman splitting B (given by B — 
gfisB/2, where B is an in-plane magnetic field, [Ib is 
the Bohr magneton and g is the wire g- factor). When 
the wire is proximity coupled to an s-wave superconduc- 
tor of gap Ag, a SC pairing term Agt < Ag may be 
induced. The wire is then arranged into a Josephson de- 
vice as illustrated in Fig. [l^,. The resulting hybrid system 
SS'NS'S can be driven into a TS phase where two MBS 
are formed at the ends of each S' region, denoted by 71,2 
and 73,4, when B > B c = y/fj? s , + A|, gj [5]. The low- 
energy Majorana sector appears for energies e < A off , 
where A e g is the effective SC gap (Fig. [2]d). 

Isolated MBS 7.; are zero energy superpositions of a 
particle and a hole. A pair of Majoranas 7^2 may be 
fused into a Dirac fermion c 1 = (71 — «72)/v / 2, such 
that the operator 2c^c — 1 = 2^72 with eigenvalues - 
1 and 1 in states |0) and |1) = c+[0) defines FP gS]. 
Thus, a spatial overlap of two MBS hybridizes them into 
eigenstates of opposite FP. In a TS wire of length Lgi, 
the decay distance of the two MBS pinned at the wire 
ends is the effective coherence length £ ef f = Hvp /7rA e g. 
Their overlap will induce a splitting ~ ±A c ffe~ L s' 
Similarly, two MBS at either side of a Josephson junc- 
tion with phase difference <f> and transparency T will 
hybridize into even/odd fermion states with energies 
~ ±A cff \/Tcos</>/2, [M HO Hi]. In the setup of Fig. 
[I] four MBS (two 1 inner' 72,3 and two 'outer' 7^4) hy- 
bridize both through the Josephson junction (region N) 
and the finite length S' regions. The resulting eigen- 
states are empty and filled states of two Dirac fermions 
d\ 2 (0), constructed as two ^-dependent (orthogonal) su- 
perpositions of the two fermions c\ n = (72 + 473)/ V2 and 
c tut = (7i+*74)/v / 2~i which are themselves obtained from 
the fusion of the inner and outer MBS, respectively. We 
denote eigenenergies by E ni „ 2 and eigenstates as \n1ri2), 
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FIG. 2. (Color online) Andreev bound states (ABS) for a 
SS'NS'S junction with normal conductance G = 0.72Go (with 
Go = 2e 2 /h), an induced gap Ag/ = 218/ieV, and wire's 
length of 4fj,m (where l ao = 216nm for a InSb wire). In 
the topologically trivial phase (a), B = 0.36meV = 0.56 c , 
whereas in the non-trivial phase (b), B = l.lmeV = l.5B c . 
The dashed line denotes the wire effective gap A e g = 130/ieV 
that separates localized ABS from the quasicontinuum. 



where n\,n<z = 0,1 are the occupations of fermions d\ 
and d\. Two of them have even total parity, | J, e ) = |00), 
I fe) = 1 1 1) = <4<4|00), and the other two are odd, 

I ; c ) = |io) = 4|oo), 1 to) = |oi) = 4|oo) m\- 

Enwzifi) anti-cross at <f> — ir within same-parity sectors 
(Fig. [l]}) and, hence, the supercurrents are 27r-periodic. 

A 47r-periodic Josephson effect can, nevertheless, be 
recovered by inducing LZ transitions with a voltage bias, 
such as I 4- e ) — *• I te) (green arrow in Fig. lb) [30] . To de- 
scribe the response of a realistic biased junction, however, 
an extension of the simplified Majorana model above is 
required. Indeed, non-adiabatic driving may also induce 
inelastic transitions into delocalized states above the TS 
gap (yellow arrow in Fig. [lj)). These latter transitions 
induce an effective parity-mixing rate (yellow wiggly ar- 
row) which couples even and odd sectors (quasiparticle 
poisoning). A proper description of such dynamics in- 
volves a calculation of all the Andreev levels (both below 
and above A e ff) coupled to the continuum (above Ag) of 
the junction, which we develop in what follows. 

Andreev levels.- The full spectrum of single-particle 
eigenstates may be obtained by diagonalizing the BdG 
equations for the geometry in Fig. la [51] . We ob- 
tain H B dG = |En(44-dn4) £n, with eigenener- 
gies £„(</>) plotted in Fig. [2^i,b for a representative junc- 
tion. The corresponding single particle excitations, called 
Andreev bound states (ABS), are defined as \n(<f))) = 
dTj(0)|fi(^>)), with \Sl((p)) denoting the ground state. In 
the short junction limit, Ln <C £ e ff, only two ABS he 
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below the effective gap A c g. Upon crossing into the 
TS phase, these two levels, initially localized around the 
N region in the non-topological phase (purple and or- 
ange curves in Fig. 2i), reconnect into the two d\^(<p) 
Majorana branches, xhe lowest Majorana branch (blue 
curve in Fig. [^Jd, level d\) has a weak cf) energy depen- 
dence, and is dominated by the outer MBS (c\ ut fermion, 



A cff e" 



), whereas the one at higher energy 



(red curve, d\) results mostly from the fusion of the in- 
ner MBS (cj„ fermion, e 2 « A off v^cos cj)/2). 

Due to the finite length Lg; > £s>, a dense set (qua- 
sicontinuum) of ABS delocalized across the TS wire ap- 
pears within a large energy window A 5 > e > A e g above 
the low-energy sector. This quasicontinuum (which ap- 
proaches a continuum as Ls> — > 00), separates the Majo- 
rana branches from the true continuum of the problem at 
e > As- Importantly, for finite transparency junctions, 
the Majorana levels detach from the quasicontinuum, as 
opposed to the topologically trivial case or a standard 
SNS junction [35] where the ABS touch the continuum 
at <j> = 0, 2ir. We have verified that the detachment, de- 
noted by 8q in Fig. [T]d, increases as the transmission T 
of the normal region becomes smaller, Sq / A c ff « 1 — \/T, 
in agreement with simpler models [El [17l [26] . 

ac Josephson effect.- We now study the Josephson 
current across the junction when biased with a voltage 
V. The bias makes the phase difference <p time depen- 
dent, <f>(t) = 2eVt/h = ujjt, where ojj is the Josephson 
frequency. We calculate the Josephson current through 
the biased junction in terms of a reduced density ma- 
trix in the Majorana sector. Notably its evolution con- 
tains parity-mixing terms as a result of the Andreev lev- 
els above A e ff that have been traced out [31^- The goal 
is to identify smoking-gun features in /(t), in the form 
of fractional frequency components, that may unambigu- 
ously determine the existence of MBS in the TS phase. 
It has been shown, however, that in the presence of fi- 
nite dissipation (e.g. non-zero parity mixing), the cur- 
rent becomes, in the long time limit, strictly 27r-periodic, 
as a consequence of Floquet theorem, and that as a re- 
sult nothing remains in the stationary current to qual- 
itatively distinguish the trivial from the TS phase [2"6] . 
The two phases, however, exhibit crucial differences in 
their spectra that give rise to very different features in 
their transient regimes, I(t < tx)- The spectral features 
of the TS phase (sharp <f> = n anticrossing and detach- 
ment of the Majorana sector from the quasicontinuum in 
the TS phase, see Fig. [2]) allows the even and odd Ma- 
jorana subspaces in the TS phase to evolve almost co- 
herently performing 47r-periodic population oscillations 
during a potentially long transient time tx 3> T 27r , where 
T 27T = 2tt/ujj is the period of the Josephson driving. We 
now characterize the transient regime in detail. At bias 
2eV < 5^, the junction will remain in the initial ground 
state 1 00), resulting in an adiabatic 27r-periodic current 
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FIG. 3. (Color online) Top: diagrams of the different dynam- 
ical regimes of the ac Josephson current as a function of junc- 
tion bias V and dimensionless conductance G/Gq, for (a) the 
topologically trivial and (b) non-trivial phases of a 8.7/xm long 
InSb wire similar to that of Fig. [2] Typical current profiles 
are presented in (c). Conventional adiabatic (blue, top) and 
multiple Andreev reflection (MAR) (green, bottom) regimes 
are both 27r-periodic. A wide region of beating 47r-periodic ac 
Josephson current (red, middle) appears in the non-adiabatic 
TS phase, and is absent in the trivial case. Each regime has 
a distinct frequency transform (d). The inset shows the beat 
splitting ±5ui of the anomalous peak at loj/2 as a function 
of the bias V, both for the even (solid) and the odd sectors 
(dashed). Panel (e) presents a typical dissipative decay from 
the 4-7r-regime into the stationary MAR regime, which devel- 
ops a dc current component (peak at uj — 0). 



(blue curve in Fig. [3J:). Increasing 2el^ above 8^ in- 
duces LZ transitions at (f> = n from | J,e) = 1 00) into 
I te) = 1 11) 1 which tend to produce a perfect population 
inversion as 2eV 3> 5-n ■ Then, by virtue of the finite de- 
tachment Sq, if U remains smaller than a certain Vmar 
(quantified below), the population escape into the qua- 
sicontinuum as <p crosses 27r will be negligible, and the 
junction will subsequently re-invert into its ground state 
I 4, e ) at 4> = 37T. The repetition of this inversion process 
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FIG. 4. Transient time scales as tr = ,H F,1-> r, with dimen- 
sionless function t(G/Go) plotted in (a). The transient 47T 
regime (red shaded region, b) becomes longer as the relax- 
ation rates to the quasicontinuum Fq increase (Quantum Zeno 
effect). The two main panels in (b) show the ac Josephson 
current vs. time (2eV = 0.03Ag/) and two values of To high- 
lighted in the inset, which shows the transient time tr oc Tq. 

gives rise to a transient 47r-periodic Josephson current 
/ ~ sin 0/2 (red curve in Fig. |3fc), with a superimposed 
beating envelope (due to imperfect LZ transitions at 7T, 
37r, etc.). The parameter window required by this solu- 
tion, S n <C 2eV <C 2cVmar, is experimentally relevant, 
since and Vmar may be independently controlled by 
the length Lg> and e.g. the junction transparency, re- 
spectively. In contrast, in the topologically trivial phase, 
only the lowest of the two Andreev levels is detached 
from the quasicontinuum. Thus, the state |11) will im- 
mediately decay into 1 10} upon crossing <fr = 2tt. As 
a consequence, the 47r-periodic current cannot develop. 
The anomalous 47r-periodic component in the transient 
regime, therefore, is an unequivocal signature of the ex- 
istence of MBS in the junction. At high V ~ Vmarj the 
transition into the continuum becomes significant in each 
period T 27r , for both phases. For such voltages, tx ~ T 27T , 
and close to two fermions escape into the contacts per cy- 
cle. This process yields a finite dc current (green curve in 
Fig. |3j;), which is the analogue of the multiple Andreev 
reflection (MAR) mechanism of conventional junctions. 
It should be noted that the current for any finite bias V 
eventually develops, in the stationary regime, a finite dc 
component, that has the same origin as MAR, i.e., the 
promotion of Cooper pairs into the quasicontinuum by 
iterated scattering processes (see Fig. [3)3). 

A parametric diagram of the MAR (green), the 27r- 
adiabatic (blue), and the 47r-beating (red) regimes pos- 
sible are presented in Fig. [3^i,b. The boundaries be- 
tween regimes are crossovers, denned by V — Vmar and 
2eV — (5^. The different regimes may be clearly dis- 
tinguished using a finite time spectral analysis of the 
current, shown in Fig. [3jd. The adiabatic regime ex- 
hibits components at u> = nujj, for integer n > 0, while 
the MAR regime develops an additional lu — peak. 



The 47r-periodic transient regime in the TS phase, on 
the other hand, exhibits a distinct half-integer spectrum 
u = (n + h)vj, together with a peak splitting due to 
the accompanying beating envelope. The splitting ±5lo 
is plotted in Fig. [3ji (inset) as a function of the applied 
bias, which reveals oscillations due to the interference of 
successive LZ processes at odd multiples of cf) = n. 

Transient times.- The characteristic decay time of 
the transient current I(t) ~ e~*/* T is given by tx = 
(2eV) 2 T ' w hi cn scales with bias V and the characteristic 
relaxation rate To that models the escape of quasipar- 
ticlcs from the quasicontinuum levels into the reservoirs 
|31j . The dimensionless function t(G/Gq), that depends 
solely on the junction's normal conductance, is numeri- 
cally computed in Fig. [4ji. It is maximum in the tunnel- 
ing regime and it decays to zero as transparency goes to 
one. This means that tx decreases when the detachment 
So —> 0, since then, the escape probability into the quasi- 
continuum in each cycle goes to one. The characteristic 
escape voltage Vmar is determined from the condition 
tx = 22^, 2cVmar = hT$T /2t:. One can see that, while 
a sufficiently slow driving V will suppress quasiparticle 
poisoning and will thus increase the duration of the tran- 
sient regime, a complementary, and potentially preferable 
route is to engineer the environment's T in order to in- 
crease the lifetime of the Majorana qubits without slow- 
ing operation. Indeed, by increasing the escape rates, 
parity mixing will be suppressed, increasing tx- This is 
known as the quantum Zeno effect, whereby a quantum 
system coupled to a bath may develop longer coherence 
times when increasing the speed at which information 
is lost into the bath. In a very fast environment, the 
transitions into the quasicontinuum at <j) — 2nn will be 
suppressed by this same Zeno mechanism. This is demon- 
strated in the simulation of Fig. [4]d, which shows how tx 
is tripled as a result of increasing the r . Using realistic 
parameters for InSb nanowires, we estimate that typical 
tx reach into the /j,s range at /zV bias voltages. 

Owing to these long transients, the spectrum of mi- 
crowave radiation should show clear features of the frac- 
tional frequencies |31j . Such measurement can be per- 
formed with an on-chip detector by using, for example, 
the photon assisted tunneling current of quasiparticlcs 
across a superconductor-insulator-superconductor junc- 
tion capacitively coupled to the TS one. Importantly, it 
has been already demonstrated (331 that such technique 
allows a direct detection of fractional Josephson frequen- 
cies which paves the way for MBS on-chip detection [31j . 
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Kouwenhoven and S. Kohler for useful discussions. We 
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ANALOGY WITH THE JOSEPHSON EFFECT 
THROUGH A SINGLE RESONANT LEVEL 
MODEL 

The bound states of a junction formed by a resonant 
level, of width T and energy e, sandwiched between two 
superconductors, are the solutions of [1] 



[to — e - 



Tto 



VA 2 - uj 2 



W + £ + 



-r 2 cos 2 



VA 2 - oj 2 
-0, (1) 



2 (A 2 - uj 2 j 
which in the resonant case, e = 0, reduces to 



bJ ± A cos ■ 



LO 



Va 2 ^ 



to 



= 0. 



(2) 



The corresponding eigenstates are not single fermions 
but rather electron- hole superpositions, namely Andreev 
bound states. The solutions of Eq. (2) become particu- 
larly simple in the limit r 3> A where the Andreev levels 
can be written as 



c^±A[l-2(^) 2 ]cos^. 



(3) 



For £ 7^ 0, the Andreev bound states can still be written 
analytically in the limit 4 — > as: 



with transmission r 



±Aa/1 



i+(f) 2 



rsm -, 



(4) 



In the general case, the Andreev levels of this problem 
exhibit anti-crossings at <j> = it for r ^ 1 (namely £ 7^ 0) 
and a finite detachment from the gap edges at <f> = 2im as 
Y changes. Interestingly, as one approaches £ — > 0, the 
junction becomes transparent r — > 1 and one gets An- 
dreev levels that cannot be distinguished from the ones 
of a Majorana junction like the one described in the main 
text. The deep difference between both systems appears 
when one takes into account fermion parity conservation 
when calculating the Josephson current. In the trivial 
resonant level model, the crossing at <p — n is not pro- 
tected. Since the bound states have not definite fermion 
number, transitions between positive and negative energy 
solutions are allowed. Thus the system evolves across 
<f> = it by remaining in the ground state. As a result, 
the Josephson current is 27r-periodic. In the Majorana 



case, the zero-energy bound state is a true fermion res- 
onance formed by fusing two Majoranas (each of them 
non-locally entangled with the corresponding one at the 
nanowire end) such that the constraint of fermion par- 
ity conservation leads to a protected crossing at cj> = tt 
and hence the system remains in an excite state after 
the crossing. This leads to a 47r-periodic current in the 
ideal infinite wire case, as discussed in the main text. We 
emphasize again that this deep difference between both 
systems originates from parity conservation and hence 
cannot be established by level spectroscopy alone. 

Similar considerations apply to junctions in which a 
Kondo resonance lies between two superconductors. In 
such case, the Kondo resonance is always pinned at £ = 
and, again, the Andreev level spectrum shows both 
crossings at <j> = ^ an d detachment from the gap [2, 3]. 



FULL ANDREEV LEVEL SPECTRUM FROM 
BOGOLIUBOV DE GENNES EQUATIONS 

The computation of the Andreev levels in an SS'NS'S 
junction is performed within the Nambu formulation of 
the junction Hamiltonian 



(5) 



-A* a ,(x) -H^,(x)* 



i><T'{x) 



where matrix H^°\x) models the one-dimensional semi- 
conducting wire plus left and right leads, without su- 
perconducting pairing. It includes a constant transverse 
spin-orbit coupling a, a Zeeman field along the wire 
B = gfisB/2 (where B is an external magnetic field, [Ib 
is the Bohr magneton and g is the nanowire g- factor), 
and a position-dependent band shift n(x) that accounts 
for a larger electronic density in the leads than in the 
wire [4, 5], 



2m* 



+ iaa y d x + Ba x - \i(x). 



We approximate fj,(x) by a piecewise constant function, 
[is* in the wire, fi^ in the central normal region (control- 
lable via a gate voltage V g ) and us in the leads, see Fig. 
la in the main text. The superconducting pairing also 



2 



varies with x, 

A aa >{x) = -iA(x)a y . 

A(x) is assumed weaker in the superconducting part of 
the wire than in the leads \As>\ < \Ag\, and zero in the 
central N region. We also denote the phase difference 
across the N region by </>, so that arg [A (±|a;|)] = ±0/2. 

To compute the Andreev spectrum as a function of <j> 
one may discretize the integral in Eq. (5), which ap- 
proximates the continuum problem by a Nambu tight- 
binding chain. Exact diagonalization of the correspond- 
ing H matrix yields H = \ J2 n {d+d n - d n d+) e n , where 
the single particle spectrum £„(</>) are the energies of 
the Andreev states \n{4>)) = d+(0)|fi(</>)) discussed in 
the main text, and \fl ((f))) is the ground state. Alterna- 
tively, one may employ a wavematching method on the 
continuum model, and find the energies that yield nor- 
malizable solutions in the absence of incoming modes in 
the leads. We have checked that both methods yield 
the same results for a sufficiently fine tight-binding dis- 
cretization, and exhibit band reconnection and Majo- 
rana branches if B > J A|, + y? s , , see Fig. 2b in the 
main text. The effective gap in the topological phase, 
A c ff, is the smallest of the two gaps Ai = \B — B c \ and 

A 2 = As I _ 2(l+A + v/r+2A+B') +0(A 2) with tilde 

quantities denoting energies in units of the SO energy, 
E so = h 2 /ml 2 , with I so = h 2 /m*a the SO length given 
in terms of the wire's effective mass m* and SO coupling 
a. 

NON-ADIABATIC DYNAMICS 

As a result of a bias V across the junction, the phase 
<f> becomes time dependent, <j>(t) = 2eVt/h = cjjt, and 
the system is driven out of equilibrium. For a given 
phase </>, one may build a set of many body instanta- 
neous eigenstates {\y> n{<!>)) = d +nN ■ ■ ■ d +ni |0)}, where 
n m = 0, 1 is the occupation of the m-th Andreev level 
and H(<f>)\<pfi(<l>)) = E R (4>)\ip n (4>)) . The many body in- 
stantaneous energy is E n (<j>) = J2 m £ m(4>)(n m - 1/2). 

The time evolution of the system's full density matrix 
p(t) may be expressed in the {|<Pn(i))} basis, Pnn'(t) = 
((ffi(t)\p(t)\(pfi'(t)), and is governed by the equation of 
motion, 

d tP (t) = - l -[H(t),p(t)], 

The effective Hamiltonian H(t) includes the connec- 
tion between instantaneous eigenstates, Aft'n{<P) = 
i(<Pn'(<l>)\d<i><Pn(<l>)), 

Unn>(t) = 6n,H>E n ((f>(t)) - H [d t <t>{t)] Ann' (4>(t)) . 

The many body connection Ann' may be related to 
the single particle connection A nm = i{n\d<f,\m) = 



i(O|d„(90d+ )\Q) [6] and the anomalous connection 
A nm = i(Q\(d^d+)d+\n) by 

i{<Pa\d<t>\<Pa) = X] ( 6 ) 

m 

i(fn\d4,dfdj\ipfi) = -i(<pn\d<i>didj\<pft)* = A 3l (7) 
i{tp n -\d<pd}d1 j \Lp i i) = -iifnld^didjlipn)* = A 3i (8) 

for any state \<Prt). All other connections are zero, which 
implies that the system evolution preserves fermion par- 
ity (i.e. if the initial number of excitations is definitely 
even or odd, it will remain that way, even though the to- 
tal number may change due to the anomalous connection 
A.) 

Driving may impart the excitations with energy, which 
may coherently promote them into the quasi-continuum 
above A e ff. To account for such quasiparticle poisoning, 
one could pursue a brute-force approach by including an 
extended basis with all many-body states generated by 
filling the first N Andreev levels in the spectrum. [7] 
This is impractical, however, since the size of the basis 
grows exponentially as 2 N . Moreover, it soon becomes 
clear that once a fermion fully escapes from the Majorana 
sector, it has a vanishing probability of returning, given 
the large phase space available in the quasi-continuum. 
Thus, these excitations are irreversible to all practical ef- 
fects. This single- fermion escape process changes parity. 
To model such two-step quasiparticle poisoning process 
without having to track the dynamics of the (exponen- 
tially) large number of many body states, one may sub- 
stitute the single-particle quasi-continuum (h, or 'high') 
levels \m)h by dissipative levels (with zero mutual con- 
nection) that decay directly into a fermion reservoir at 

a rate while the 'low' levels (I, below the gap) re- 

main non-dissipative. If we assume a Markovian approx- 
imation for the decay of h into the fermionic bath, the 
master equation takes on a Lindblad form 

dtp = -\ [H, p]+J2 U mP L+ \ {L+L m ,p}) 

If we further constrain our dynamical space to states 
with a total of I or fcrmions in the h sector, 
{|<Pn ( ;i m )) |Vni;o)} (fast decay limit), the Lindblad oper- 
ators L m will simply project any state with a fermion in 
the m-th dissipative level (m € h), into another without 

it L m = J2n t l^« ( :0)(^nt;l m |- 

For large enough <~ O (v 1 ) (v being a small per- 
turbativc parameter), the irreversible decay will suppress 
the density matrix in the quasi-continuum, pu <~ 0(1), 
plh,Phl ~ 0{r}), phh ~ O {rj 2 ). One may then pertur- 
batively solve the detailed balance conditions for the h 
sector, which eventually leads to a Lindblad-type master 
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equation for the reduced density matrix p a , 
dtPu = [Wu,Pu] 

+ ^r a/3 (c aPll c+- l -[c+c a:Pu ])+o(T 1 ) 

a/3 ^ ' 

Considering the I sector as spanned by two Andreev lev- 
els, of energies £1,2, the corresponding many body ba- 
sis |<Pni,n 2 ;0) = l"-1^2> IS {| 4-e), I te), I io},\ to)} = 

{|00), |11), 1 10), |01)}, where e and o stand for even and 
odd fermion parity. Then, the Lindblad equation above 
reduces to 



dt~P~-\ [H,p[ +J2 T <*0 UoPtf - I {C+C a ,p}\9) 

a/3 ^ ' 

with four parity-mixing Lindblad operators, 



A= Ue)(to| + Uo)(te 
C 2 = I ie)(io I - I to)(te 
r 3 = -|te)(to| + Uo)Ue 
A = |te)Uo| + |to)Ue 



(10) 



and a 4 x 4 relaxation matrix 



r a/3 («/>) = 4c4£ 



{4>)v m p{4>) 

r (m) 



(11) 



The v m {<j)), given explicitly in terms of single-particle 
connections between the I states and the h states, 



V-n — {j42.ru Ai. n , Ai^ n , j42. n } , 



(12) 



are peaked at around <f> = 2im, for integer n. They 
quantitatively account for parity mixing mediated by 
coherent excitation into dissipative level m, and con- 
tain detailed microscopic information about the quasi- 
continuum, making the resulting dynamics for the re- 
duced density matrix highly non-trivial. 

Since the levels e m >3 above A e ff are almost <f> inde- 
pendent, the Josephson current through the biased junc- 
tion may be approximated in terms of the 4 x 4 re- 
duced density matrix in the Majorana sector, namely the 
P„ 1 n 2 ,n' 1 ni = { n i n 2 |pK™ 2 ) governed by Eq. (9). Then 



I(t) = 



4e x - 

«1,«2 = {0,1} 



Pn 



i«2 ,nm2 



(*) fy£ nina (0(t)),(13) 



where the many body energies E nin2 are the Majorana 
branches of Fig. lb in the main text. Solving p(t) with 
Eq. (9), we compute the Josephson current under bias 
V, assuming the junction is initially in its ground state 
|ft(0)) (at time t = 0). 



DISCUSSION ABOUT EXPERIMENTAL 
DETECTION 

Parity protection by the Zeno effect should prove use- 
ful for prolonging the 47r-periodic transient regime when 
the main source of quasiparticle poisoning is through the 
quasi-continuum. Other sources of poisoning not dis- 
cussed here include non-equilibrium quasiparticles from 
the rest of the circuit, a contrition that can, however, 
be greatly reduced by including quasiparticle traps, such 
as nearby metal contacts [8] . In the case of quasicontin- 
uum poisoning, a rough estimate relating the quasiparti- 
cle escape velocity to the Fermi velocity yields r « A e g, 
which suggests that increasing A e ff will improve Zcno- 
type parity protection of driven Majorana qubits. Using 
realistic parameters for InSb nanowires, we estimate that 
typical transient times reach into the /xs range at pN bias 
voltages. 

The spectrum of microwave radiation from TS junc- 
tions should show clear features of the fractional frequen- 
cies (Fig. 3d, main text). Such measurement can be 
performed with an on-chip detector, which greatly mini- 
mizes the impedance-matching problems in classical de- 
tection schemes. This on-chip detection can be achieved 
by using, for example, the photon assisted tunneling cur- 
rent of quasiparticles across a superconductor-insulator- 
superconductor junction capacitively coupled to the TS 
one. It has been already demonstrated [9] that such 
technique allows a direct detection of fractional Joseph- 
son frequencies [10] of a superconducting single electron 
transistor [11]. This is possible if 7\ 3> tx, where T\ 
is the relaxation time corresponding to the transition 
I t e ) = 1 11) — > I l e ) = 1 00) that brings the system back 
to the ground state. Such parity-conserving transitions, 
that at low enough temperatures are due to quantum 
fluctuations of the phase, can be attributed to photon 
emission to the electromagnetic environment and, to a 
lesser amount, to phonon emission. Proper engineering of 
the circuit containing the topological wire can, in princi- 
ple, greatly reduce both emission processes. For example, 
photon emission can be inhibited by reducing the ohmic 
component of the impedance seen by the Majorana qubit 
[12]. To lowest order, this relaxation rate reads: 

■ ' | ( t.|iu.)| ^' z ' AE ^" ,») 



AE e 4e 2 
8AE, 



R 



Kte |^Ue)| 



K 

2 Re[Z(AE e /h)} 



R 



K 



where I = (4e/h)d^H, AE e = E te - E le and R K = 
h/e 2 w 25.81fcf2 is the resistance quantum. This rate 
is (/(-dependent. Substituting the value of the connec- 
tion between the even Majorana branches and their en- 
ergy difference for the wire considered in Fig. 3 of the 
main text, and averaging over the phase, we obtain re- 
laxation times of the order of Ips for impedances of 
Re [Z(AE e /h)\ = m. 
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